Quantum phase transition in the one-dimensional compass model 
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We introduce a one- dimensional model which interpolates between the Ising model and the quan- 
tum compass model with frustrated pseudospin interactions a^af^i and afaf^i, alternating between 
even/odd bonds, and present its exact solution by mapping to quantum Ising models. We show that 
the nearest neighbor pseudospin correlations change discontinuosly and indicate divergent correla- 
tion length at the first order quantum phase transition. At this transition one finds the disordered 
ground state of the compass model with high degeneracy 2 x 2^''^ in the limit of — > oo. 
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I. INTRODUCTION 

Recent interest in quantum models of magnetism with 
exotic interactions is motivated by rather complex su- 
perexchange models which arise in Mott insulators with 
orbital degrees of freedom.^ The degeneracy of 3d orbitals 
is only partly lifted in an octahedral environment in tran- 
sition metal oxides, and the remaining orbital degrees 
of freedom are frequently described as T — 1/2 pseu- 
dospins. They occur on equal footing with spins in spin- 
orbital superexchange models and their dynamics may 
lead to enhanced quantum fluctuations near quantum 
phase transitions,^ and to entangled spin-orbital ground 
states.- The properties of such models are still poorly 
understood, so it is of great interest to investigate first 
the consequences of frustrated interactions in the orbital 
sector alone. 

The orbital interactions are intrinsically frustrated — 
they have much lower symmetry than the usual SU(2) 
symmetry of spin interactions, and their form depends 
on the orientation of the bond in real space,— so they 
may lead to orbital liquid in three dimensionsi^ Although 
such interactions are in reality more complicatedj^i^i^i^ a 
generic and simplest model of this type is the so-called 
compass model introduced long agOf^ when the coupling 
along a given bond is Ising-like, but different spin com- 
ponents are active along different bond directions, for in- 
stance JxO'fo'j and JzCrfcr^j along a and h axis in the two- 
dimensional (2D) compass model. The compass model 
is challenging already for the classical interactions.^ Al- 
though the compass model originates from the orbital 
superexchange, is is also dual to recently studied models 
oi p -\- ip superconducting arraysi^ It was also proposed 
as a realistic model to generate protected cubits,^ so it 
could play a role in the quantum information theory. 

So far, the nature of the ground state and pseudospin 



correlations in the compass model were studied only by 
numerical methods. It has been argued that the eigen- 
states of the 2D compass model are twofold degeneratci^ 
In contrast, a numerical study of the 2D compass model 
suggests that the ground state is highly degenerate and 
disordered 3 In fact, the numerical evidence suggests a 
first order quantum phase transition at Jx = Jz-^^ with 
diverging spin fluctuations and a discontinuous change in 
the correlation functions. 

The purpose of this paper is to show by an exact solu- 
tion a mechanism of a first order quantum phase transi- 
tion in quantum magnetism. Such a transition from qua- 
siclassical states with short range order to a disordered 
state occurs in the one-dimensional (ID) XX-ZZ model, 
with antiferromagnetic interactions'^^ between erf and al 
pseudospin components, alternating on even/odd bonds. 
The model is solved analytically by mapping its different 
subspaces to the exactly solvable quantum Ising model 
(QIM)^ which plays a prominent role in understanding 
the paradigm of a quantum phase transition, including 
recent rigorous insights into the transition through its 
quantum critical point at a finite rateji^ Note that the 
model introduced below is generic and by no means lim- 
ited to the orbital physics. For instance, the QIM is 
realized by the charge degrees of freedom in NaV2 05, 
where it helped to explain the temperature dependence 
of optical spectra.iS. 

The paper is organized as follows: In Sec. II we intro- 
duce the pseudospin (orbital) model which realizes a con- 
tinuous interpolation between the classical Ising model 
and the compass model in one dimension. This model 
is next solved exactly in its subspaces which are sepa- 
rated from each other. The properties of the model and 
the mechanism of the quantum phase transition are eluci- 
dated in Sec. III. At the transition point we demonstrate 
the discontinuous change of correlation functions (Sec. 



2 



III. A) and the vanishing pseudospin gap (Sec. III.B). 
The paper is summarized in Sec. IV. 



II. PSEUDOSPIN XX~ZZ MODEL 

In order to understand the nature of a quantum ground 
state with high degeneracy found in the ID compass 
model with alternating superexchange interactions be- 
tween erf and erf pseudospin components, it is important 
to investigate how the pseudospin correlations develop 
when this state is approached. Therefore, we start below 
from the classical Ising model with interacting erf pseu- 
dospins and modify gradually the interactions on every 
second bond by replacing them by the ones between erf. 
Although this does not correspond to any deformation of 
interacting orbitals in a crystal, in this way we keep a 
constant value of the total pseudospin coupling constant 
J on each bond, and can demonstrate that the frustra- 
tion of interactions increases and dominates the behavior 
of the compass model at the quantum critical point. 



A. Hamiltonian and its subspaces 

We consider the pseudospin model with interactions 
depending on parameter a S [0, 2]. When < a < 1 

JV' 
i=l 

(1) 

with N' = N/2, but when l<a<2 

N' 

+ (2-aK.<+i + (a-lK,a2Vi]. (2) 

Here {erf ^al} are Pauli matrices. We assume even num- 
ber of pseudospins N and periodic boundary conditions: 
cn+1 = o"!. The model reduces to the ZZ Ising model at 
a = 0, 

JV 

H(0) = J^afaf+i, (3) 

i=l 

and to the XX Ising model at o; = 2, 

N 

W(2) = J^«+i. (4) 

i=i 

It is frustrated between these two incompatible types of 
order for < o; < 2, and we focus on the consequences of 
increasing frustration. Right in the middle between the 
two classical cases ^ and i.e. for a = 1, it becomes 



the proper quantum Y^-TIL model (called ID compass 
model), 

N' 

7i(l) = jJ^hViaf. + a^^Vi] . (5) 

favoring antiferromagnetic order of x and z pseudospin 
components on odd and even bonds, respectively. 

The Hamiltonians H]) and ^ are related by symmetry: 
simultaneous exchange of erf erf , 2i <-> 2i — 1 V i, and 
(1 — a) <-*■ (of — 1), maps the Hamiltonians into each 
other. Thanks to this symmetry, we need to solve only 
the model ([T]) for o; € [0, 1], with modulated interactions 
for odd pairs of pseudospins {2i — 1, 2i} (on odd bonds). 
The Hamiltonian ([T]) can be conveniently solved in the 
basis of eigenstates of the ZZ Ising model at o; = 0. We 
start with one of the two degenerate ground states, 

inTJ.TJ.TJ.TJ.Ti...) , (6) 

which is modified at finite a > by spin-flip operators: 
— rf = o-2i_ierfj. Action of various r^'s on the state © 
generates other eigenstates: 

liTTiTiTiTiTi...) , 

ITiiTTiTiTiTi...) , ... , 

liTiTTiTiTiTi...) , ... , 

liTiTiTiTiTiT...) , (7) 

where the last state is the second degenerate ground state 
of the ZZ Ising model. This set of states (H))-© is a 
convenient basis in the subspace where all odd pairs of 
pseudospins are antiparallel, and give a constant energy 
contribution of —(1 — a)J per bond. By construction, 
the Hamiltonian ([1]) does not mix this subspace with the 
rest of the Hilbert space. Therefore, in this subspace the 
energy due to the even {2i, 2i + 1} bonds in Eq. ([T]) 
depends on the pseudospin orientation on the two neigh- 
boring odd bonds, and may be expressed by a product 
T^T^j^i of two spin operators, with — rf = crfi-icl;. 
this representation it is clear that the Hamiltonian ^ 
reduces in this subspace to the exactly solvable QIM.^'^ 
In general the Hilbert space can be divided into sub- 
spaces where different odd pairs of spins are either par- 
allel or antiparallel. Each subspace can be labelled by a 
vector s = (si, sni), with Si = 1 (si = 0) when two 
pseudospins of the odd bond {2i — 1, 2i} are parallel (an- 
tiparallel). In each subpace s the terms oc (1 — a) for odd 
bonds in Eq. ([T]) give a constant energy contribution 

N' 

Cs{a) = {l-a)JY,^2^-1^2^ 

1=1 

= -{l-a)J{N' -2s), (8) 

where s = X^i^i number of parallel odd pairs of 

spins. With a convenient choice of spin operators for the 
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antiparallel {2i — 1, 2i} pairs of pseudospins, 

-Tf = (lTi)(iTl + liT)(ni) , 

-rf = (inxni-UTXiti) , (9) 

and for the parallel pairs 

= (ITT)ail + lii)(TTI) , 

-rf , (lo) 

the Hamiltonian ([T]) reduces to 

1=1 

- J [«ri^, + (-1)^ T^,Tf ] + C,(a) . (11) 

For even s the Hamiltonian (jlip is simply the ferromag- 
netic QIM, but when s is odd, the interaction on the last 
{N', 1} bond is antiferromagnetic. 

B. Exact solution 

Each effective model (fTTj) can be solved using the 
Jordan- Wigner transformation for spin operators, 



(c. + 4)n(i-24s 



(12) 
(13) 



Here Ci is a fermionic annihilation operator at site i. Af- 
ter this transformation the Hamiltonian pip becomes 



H,ia) = J E ( 

i=l 



ac]c 



c a 



i H+l 



+ h.c.) 

h.C.) 



with 



Cl 



^ (-1)1+^+^^ 



(14) 



(15) 



are projectors on the subspaces with even 
(— ) numbers of c-fermions and 



and odd 



+ Cs{a) (18) 

are corresponding reduced quadratic Hamiltonians. Def- 
inition of the Hamiltonians (jlSp is not complete without 
boundary conditions which depend both on the choice 
of subspace ± and on the parity of s, see Eq. (|15l) . 
When both s and the number of c-fermions have the same 
parity, then the boundary condition is antiperiodic, i.e., 
cat'+i = — Cl, but when on the contrary the two numbers 
have opposite parity — it is periodic, i.e., ctv+i = C\. 
With these boundary conditions the Hamiltonians (fTS]) 
and (|14p are the same. 

The Hamiltonian (|18p is simplified by a Fourier trans- 
formation, Cj = X]fc '^fc^'*^"' • Here fc's are quantized 
pseudomomenta. In the following we assume for con- 
venience that N' is even (i.e., N = 4m, m integer). 
For periodic boundary conditions fc's take "integer" val- 
ues (recall that N = 2N') k = 0,±j^,±2j^, ■ ■ ■ , tt , 
and in the antiperiodic case they are "half-integer", i.e., 
k = ±i|f^,±f^,-- - ,±^{N' . As a result, 

the Hamiltonians describe independent subspaces 
labelled by k, with mixed k and — fc quasiparticle states. 

Hp (a) — ^2{a - cos fc)c^c^ -I- sin A; (^cj^clj, + h. 



+ Csia) 



(19) 



Diagonalization of Hi{a) is completed by a Bogoli- 

ubov transformation,"'^'^ Ck = Ukjk + V-k^-k' where 
the Bogoliubov modes {uk,Vk} are eigenmodes of the 
Bogoliubov-de Gennes equations: 

euk = 2J{a — cosk) Uk + 2J sink Vk , (20) 
evk = 2J sink Uk ~ 2J{a — cosk) Vk ■ (21) 

For each value of k there are two eigenstates with eigenen- 
ergies Cfe and — e/c- Positive eigenenergies 



depending on parity of the number of c-quasiparticles. 
This parity is a good quantum number because c- 
quasiparticles can only be created or annihilated in pairs, 
as can be seen in Eq. (fT4|) . 

In order to extend the sum in the first line of Eq. (fT4|) to 
i = N' and, at the same time, to include the second line 
into this extended sum, we split the Hamiltonian (|14p as 



H^ia) = P+ Hi (a) P+ + p- HZ (a) p- 
where 



N' 


1 

^ 2 




i=l 





N' 



2c]c, 



(16) 



(17) 



= 2J\J\ + a2 _ 2a cos /c 



(22) 



define quasiparticles 7^ for each k which bring the Hamil- 
tonian (fT9|) to the diagonal form 



(23) 



being a sum of fermionic quasiparticles. However, thanks 
to the projection operators P^ in Eq. (fT6|) only states 
with, even or odd numbers of Bogoliubov quasiparticles 
belong to the physical spectrum of Hg{a). 

In order to find if the number of Bogoliubov quasipar- 
ticles in a given subspace must be even or odd we must 
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FIG. 1: (Color online) Intersite pseudospin correlations on 
odd {2i - l,2i} and even {2i,2i + 1} bonds in the XX-ZZ 
model H}-© for increasing a. A crossover between two 
types of pseudo-order, with (cr|i_i(T|i) — —1 for a < 1 and 
{(Tfj(T|i_i_i) = —1 for a > 1, occurs at the quantum critical 
point a = 1, where only (crfi_icrfi) = {cr|icr|i+i) = -| are 
finite. 



find first the parity of the number of c-quasiparticles in 
the Bogoliubov vacuum in this subspace. This parity de- 
pends on the boundary conditions. Indeed, when they 
are antiperiodic, then sin k ^ for any allowed k and 
the ground state of the Hamiltonian p9)) is a superposi- 
tion over states with pairs of quasiparticles c^cL^j. This 
Bogoliubov vacuum has even number of c-quasiparticles. 
On the other hand, for periodic boundary conditions we 
have two special cases with sin fc = when fc = or 
IT. When < a < 1 the Bogoliubov vacuum contains the 
quasiparticle cq but not the quasiparticle c^r and the num- 
ber of c-quasiparticles is odd. In short, (anti-)periodic 
boundary conditions imply (even)odd parity of the Bo- 
goliubov vacuum. 

Taking further into account that any 7^ changes par- 
ity of the number of c-quasiparticles, we soon arrive at 
the simple conclusion that even (odd) s implies that only 
states with even (odd) numbers of Bogoliubov quasiparti- 
cles belong to the physical spectrum of Hg{a). Note that 
for odd s the lowest energy state is not the Bogoliubov 
vacuum but a state with one Bogoliubov quasiparticle of 
minimal energy t}~ . Its pseudomomentum is /c = in a -I- 
subspace, but in a — subspace there is a choice between 
two degenerate quasiparticle states with fc = ±i^. 



III. QUANTUM PHASE TRANSITION 

A. Correlation functions 

For any a G (0,1) the ground state is found in the s — 
(-I-) subspace. Pseudospin correlators in the ground state 
can be expressed by the spin correlators of the effective 
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FIG. 2: Distance dependence of {o'2i<^2{i+r)) correlator be- 
tween even pseudospins belonging to different odd bonds. 
Two limits a — > 1* demonstrate the discontinuity at the 
quantum critical point (a = 1). The correlators for a — > 1~ 
approach the asymptotic value for large r in algebraic rather 
than exponential way indicating divergent correlation length. 



QIM Eq. (fTH]) . For odd pairs of pseudospins: 



N' ^ 

fe 



(24) 
(25) 



Surprisingly, in spite of decreasing interaction oc {1 — a) 
in Eq. H]), the odd {(^ii-i'^ii) correlator ([24|) shows the 
same perfect pseudospin order as that found at a = 0, 
while the (crfi^icrfj) one (^5]) gradually decreases from 
at a = to — - when a ^ 1~ (at ^ 00), see Fig. [T] 
For different odd pseudospin pairs one finds (a < 1): 



-„) - , 

-„) = (-!)"+"< 



(26) 
(27) 



when i ^ j and for m, n = 0, 1. As is well known in the 
QIM, the correlator on the right hand side of Eq. (|27|) is 
Toeplitz determinant,— 



/l /2 
/o fl 



fr 
fr-1 



(28) 



with constant diagonals fr{a) = Sr,o — 2ar{a) + 2br{a), 
given by ar{a) = 1^ Sfe kfcP cos(fcr) and bria) = 
Ufcf^ sin(fcr). This correlator is positive for all 
r and finite when r —>■ 00, indicating long range ferro- 
magnetic order of moments. When r ^ cxd it decays 
exponentially towards finite long range limit with a cor- 
relation length ^ which diverges as ^ (1 — a)"^. As 
we have seen, in this limit the gap in the quasiparticle 
energy spectrum ([^ tends to and a quantum critical 
point is approached in the universality class of the QIM. 

For nearest neighbor even pairs [when j — i -I- 1, 
m — and n = 1 in Eq. (HZ])], the (CT|jCr|j_|_i) pseu- 
dospin correlator is negative, see Eq. ([?7|) . as expected 



5 




FIG. 3: Eigenenergies E„ of the XX-ZZ model (H])-© as 
obtained for a chain of A'^ = 8 sites with periodic boundary 
condition for increasing a. Level crossing at a = 1 marks the 
quantum critical point of the compass model ^ . 



for an antiferromagnet. This correlator shows a com- 
plementary behavior to that of the (cfi-i'^'fi) correla- 
tor for odd bonds ([25]) — it interpolates between —1 at 
(2 = and when a ^ 1~ (see Fig. [T]). In con- 
trast, the (cr|j_icr|j) and (CTfjCrfi+i) correlators change 
in a discontinuous way at the quantum critical point 
(a ~ 1), where the pseudospins become entirely disor- 
dered, and (cr|i_icr|,) = {(J^^a^^j^^) = 0. Therefore, only 

{'^2i-i'^2i) — i'^2i^2i+i) — finite and contribute 

to the ground state energy of the compass model. 

Using the symmetry of the model ((T])~(l2|), the ground 
state correlators for a € (1,2] can be obtained by map- 
ping the correlators for a € [0, 1). In this way we obtain 
correlators that are well defined for any value of a except 
a = I, where most of the correlators are discontinuous 
(Fig. [21). For example, when a < 1 we have antiparallel 
odd pairs of pseudospins ((o'fi-iO'li) = but the same 
correlator tends to when a 1+ (Fig. [1]). This discon- 
tinuity of the correlators at the quantum phase transition 
is a manifestation of level crossing between ground states 
in different subspaces swhen a ^ 1^. 



B. Ground state energy and excitations 

The spectrum of the Hamiltonian ([l} changes quali- 
tatively from a ladder spectrum with the width of 2JN 
at a = to a quasicontinuous spectrum with the width 
reduced to | JiV at the a = 1 point (Fig. [3]). While the 
ground state has degeneracy ci = 2 at a = (the -I- and 
— lowest energy states in the subspace with s = 0) and 
(i=lforO<Q;<l (the + ground state with .s = 0), 
large degeneracy occurs at a = 1. 




FIG. 4: (Color online) Pseudospin excitation gap A in the 
XX-ZZ model (H}-© for increasing a (solid lines). The gap 
collapses at the quantum critical point (QCP) (in the compass 
model obtained at a = 1), which separates the disordered 
phase with finite pseudospin (o"2iC"2{i+r)) correlations (left) 
from the one with finite {'72i'^2{i+r)) correlations (right). 

Lowest energies in different subspaces s are: 

k 

= £+ia) + 0{l), (29) 

where for a < 1 

£+{a) = -i(l-a)JiV 

- JN— [ dk^/l + a^ -2acosk (30) 
27r Jo 

is the energy of the ground state found in the s = (+) 
subspace in the limit of — > oo. A similar formula to 
Eq. ((30)) is easily obtained for a > 1 by a substitution 
a — > (2 — a). For any s the gap between the -t- and — 
lowest energy states is 0{1/N), but the gaps 2(1 — a)s 
between s = and s > are 0{1) and survive when 
N —^ oo. However, when a ^ 1^, then all ground state 
energies for +/— and for any s become degenerate, level 
crossing between all 2 x 2^/^ lowest energy states oc- 
curs when N oo. This ground state degeneracy is 
much higher than 2x2^ found in the L x L 2D com- 
pass model, but the overall behavior is similar — it 
indicates that frustrated interactions, the discontinuity 
in the correlation functions, and the pseudospin liquid 
disordered ground state are common features of the ID 
and 2D compass models. 

The ground state energy Sgia) (|30p increases with a 
for a G (0, 1) (Fig. [3]), as pseudospin interactions are 
gradually more frustrated when a ^ 1.— The ground 
state is separated from the lowest energy pseudospin ex- 
citation by a gap A = 2eo = 4J|1 — a|, which vanishes at 
a — 1 (see Fig. |4]). Note that this excitation corresponds 
to reversing a pseudospin component for a < 1, and 
a erf pseudospin component for a > 1. We have veri- 
fied that the linear decrease of A when a ^ 1 is well 
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reproduced by the energy spectra obtained by exact di- 
agonalization of finite systems (as obtained, for instance, 
from the data of Fig. [3]). 

IV. CONCLUSIONS 

We have shown by an exact solution of the ID XX- 
ZZ model ([T])-® that pseudospin disordered states are 
triggered already by an infinitesimal admixture of the 
interactions between other spin components than those 
used to construct classical Ising models at a = jS]) or 
a = 2 The properties of the XX-ZZ model are sum- 
marized in Fig. m with two different types of pseudospin 
correlations, dictated by the 'dominating' interactions. 
These pseudospin correlations and finite excitation gap 
may be seen as precursor of the antiferromagnetic or- 
der induced by the type of interactions in the respective 
classical Ising model, with either erf or erf pseudospins, at 
a = or a = 2, respectively. These opposite trends be- 
come frustrated in the ID compass model lying precisely 
at the quantum critical point. We anticipate that a sim- 
ilar quantum critical point determines the properties of 



the orbital liquid state in a 2D compass model.-lS 

In conclusion, the ID XX-ZZ model provides a beauti- 
ful example of a first order quantum phase transition be- 
tween two different disordered phases, with hidden order 
of pairs of pseudospins on every second bond. When the 
pseudospin interactions become balanced at the quantum 
critical point, the pseudopspin (orbital liquid) disordered 
state takes over. It is characterized by: (i) high 2 x 2^/^ 
degeneracy of ground state, and (ii) the gapless excita- 
tion spectrum. 
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